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We study the non-equilibrium dynamics of domain walls in real time for 4 and Sine Gordon 
models in 1 + 1 dimensions in the dilute regime. The equation of motion for the collective 
. coordinate is obtained by integrating out the meson excitations around the domain wall to one- 

loop order. The real-time non-equilibrium relaxation is studied analytically and numerically 
to this order. The constant friction coefficient vanishes but there is dynamical friction and 
relaxation caused by off-shell non-Markovian effects. The validity of a Markovian description is 
studied in detail. The proper Langevin equation is obtained to this order, the noise is Gaussian 
and additive but colored. We analyze the classical and hard thermal loop contributions to the 
self-energy and noise kernels and show that at temperatures larger than the meson mass the 
hard contributions are negligible and the finite temperature contribution to the dynamics is 
governed by the classical soft modes of the meson bath. The long time relaxational dynamics 
is completely dominated by classical Landau damping resulting in that the corresponding time 
scales are not set by the temperature but by the meson mass. The noise correlation function 
and the dissipative kernel obey a generalized form of the Fluctuation-Dissipation relation. 
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1 INTRODUCTION AND MOTIVATION 



Kinks and domain walls play a fundamental role in the equilibrium and non-equilibrium processes 
after phase transitions to broken symmetry states in theories with scalar order parameters. In scalar 
field theories that undergo a symmetry breaking phase transition the process of phase ordering pro- 
ceeds by the formation of domains of the ordered phase separated by domain walls. These domain 
walls are topological defects that separate regions in which the order parameter is locally ordered 
and therefore locally the system is in a broken symmetry ground state in each domain. Interest 
in the dynamics of these topological excitations is interdisciplinary. In condensed matter systems 
solitons (or kinks) are collective excitations in quasi-one dimensional charge density wave systems 
and conducting polymers [jl], |2|, [J. In particle physics domain walls in the form of sphaleronsQ have 
been argued to play an important role in baryogenesis||, and in cosmology domain walls and other 
topological defects could be responsible for structure formation [||, M. 

The classical statistical mechanics of a gas of kinks in one spatial dimension has been previously 
studied || and it was established that the kink density at a temperature T is approximately given 
by riK ~ e~ M l T with M being the kink mass. Therefore a study of the dynamics of single domain 
walls or kinks will be valid in the dilute regime M » T in which the interaction between kinks 
can be ignored because the mean separation between kinks is much larger than the typical width 
of a kink (of the order of the zero temperature correlation length or inverse meson mass). 

The focus of this article is to study the real time dynamics of relaxation of domain walls (kinks) 
in 1 + 1 dimensions via the interaction between the domain wall and the meson fluctuations in model 
field theories. In particular we study in detail scalar cp 4 and Sine-Gordon kinks in the dilute regime 
in which T « M. This problem is important in particle physics, condensed matter and cosmology. 
In particle physics dissipative processes on the dynamics of sphaleron configurations are important 
to establish corrections to the proper transition rates. In quasi-one dimensional condensed matter 
systems kinks and domain walls are responsible for important transport phenomena and therefore 
a study of the dissipative aspects will provide a deeper understanding of these phenomena. In 
cosmology the evolution of domain walls or dynamic of interfaces determines the scales in which 
ordering of horizon sized regions occur. 

Although a study in 1 + 1 dimension may not be a proper realization of the 3 + 1 dimensional 
situations in particle physics and cosmology, it will at least highlight important aspects of the 
dynamics that must be generalized to the proper situations. 

In condensed matter there is a considerable effort in understanding dissipative aspects of solitons 
starting from a microscopic description @-[[l^] in terms of Mori's formulation of linear response, and 
more recently in terms of a system-bath formulation) 13]. 



Recently Khlebnikov[14| has studied the velocity of a bubble wall in the case of a non-degenerate 
scalar potentials. The bubble-wall velocity was related to the self-energy of the scalar field through 
the fluctuation-dissipation theorem. Assuming a trilinear coupling to another massive field a local 
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friction coefficient was extracted. Alternatively, Arnold [|l5| provided an equivalent result to the 
one obtained in Ref. [14| at one loop using reflection and transmission coefficients for particles 
scattering off the bubble wall. 



Our approach is rather different. It is tailored to obtain a real time description of the dissipative 
processes and a consistent derivation of the Langevin equation in a weak coupling perturbative 
expansion. The main ingredient is the collective coordinate quantization of the kink that allows 
to obtain the non-equilibrium generating functional for the collective coordinate by integrating 
out the meson degrees of freedom, i.e. the fluctuations around the kink. The resulting Langevin 
equation allows an unambiguous identification of the dissipative kernel and the noise correlation 
function that obey a generalized fluctuation dissipation relation. The dissipative processes arise 
from the interaction between the collective coordinate and the orthogonal fluctuations around the 
kink, rather than from the coupling to other fields. 

We provide analytic and numerical study of the solutions of the equations of motion of the 
kink collective coordinate in lowest order (one loop) and establish that a Markovian approximation 
fails to describe the dynamics at large temperatures. Furthermore we analyze in detail the high 
temperature low density regime in which m « T « M with m being the meson mass focusing 
on the classical and hard thermal loop contribution to the dissipative kernel and the noise-noise 
correlation function. We argue that in lowest order in perturbation theory, the long time dynamics 
is completely dominated by classical Landau damping. 

The main results of this article are the following 

• A field-theoretical derivation of the real-time non-equilibrium equations of motion of the col- 
lective coordinate associated with translations of the domain wall and its solution in relevant 
cases for the sine-Gordon and </> 4 potentials. 

• A detailed microscopic derivation of the non-equilibrium influence functional, the quantum 
Langevin equation and the generalized fluctuation dissipation theorem to one-loop order. 

• A detailed analytic and numerical study of the relaxation in the one-loop approximation. The 
Markovian approximation is compared to the "exact" dynamics in a wide range of temperature 
and the high temperature and classical limits analyzed in detail. The long time dynamics to 
this order is analyzed both analytically and numerically. 

To our knowledge these aspects of domain wall dynamics had not been studied previously. 

Section || summarizes briefly the main concepts in collective coordinate quantization that are 
relevant for our study. In section || we introduce the main tools of non-equilibrium field theory to 
study the kink in a bath of mesons in equilibrium and describe in general the relevant interactions, 
the equation of motion of the collective coordinate a Markovian approximation and the Langevin 
equation. Section || analyzes in detail the </> 4 and Sine-Gordon models for which a Markovian 
approximation is shown to fail at large temperatures. In section [| we study the high temperature 
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but low density limit {m « T « M) and establish that the long time dynamics is dominated 
by classical Landau damping processes. In section || we discuss higher order corrections to the 
results obtained within the one-loop approximation and we comment on generalization to higher 
dimensions. Section presents our conclusions. Several appendices are included for technical 
details, in particular appendix]^ establishes the generalized fluctuation-dissipation relation between 
the damping kernel and the noise-noise correlation function. 



2 COLLECTIVE COORDINATE QUANTIZATION 

To begin our study of the dynamics of kinks we focus on 1 + 1 dimensional quantum field theories 
described by Hamiltonians of the form 

»-/*{t + 5(£)'-H « 

in which the potential U((f>) admits degenerate, broken symmetry minima. 
A static kink is a solution of the time independent field equation 



with boundary conditions such that 4> s (x — > dboo) = 4>±oo andC/(0 = |- oo ) = |lq]-[26|. Translational 
invariance implies that such solution is of the form 4> s (x — xo) with xq an arbitrary translation 
chosen such that 4> s (0) = 0, therefore xq is identified with the position of the kink. 



Lorentz invariance results in that a kink moving with constant velocity is given by (f> s 
[19]. The mass of the kink, i.e, the energy of a static kink is given by 



X— Xq — vt 
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M = E[<j> 8 ] = J dx f 1 . (3) 

Quantization around the static kink solution implies writing 

(f)(x,t) = 4> s {x - Xq) + ip(x - x ;t). (4) 

Where the fluctuation operator is expanded in terms of a complete set of harmonic modes around 
the kink 

oo 

■4>(x - x ;t) = ^2q n (t)U n (x - x Q ) (5) 
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where the mode functions U n (x — xq) obey 



dx 2 dd> 2 



U n (x - x ) = oj 2 n U n (x - x ) 



(6) 



with the completeness relation given by 

^ U^{x — xo)Ub(x' — xq) + / dk U£(x — xo)Uk(x' — xq) = S(x — x') (7) 
b 

and the subscript b stands for summation over bound states and k for scattering states. For bound 
states, the eigenvectors are chosen to be real and for scattering states, we label them as U}~(x — xq) 
and are chosen such that lit = U^ki m which case the coordinate operators obey the hermiticity 
condition q%(t) = q-k(t). 



These eigenvectors are normalized as 



J dxU*(x — xq)U, 



q (X - X ) 



J P,Q- 



(8) 



As a consequence of translational invariance, there is a mode with zero eigenvalue given by|16] 



an 



Uq(x - Xq) 



1 



M V dx 



0) 



Depending on the particular form of the potential U(cp) there may be other bound states (as is 
the case with the <fi A potential). There is a continuum of scattering states with frequencies io\ = 
k 2 + uj 2 ; u>1 = d 2 U((p)/d 2 (p\ ( j >ao . These scattering states correspond asymptotically to phase shifted 
plane waves in the cases under consideration because the relevant potentials are reflectionless| 1£ , 17] . 
The continuum states are identified with meson states, whereas bound states (other than the zero 
mode) are identified with excited states of the kinkll^ 



The fluctuation along the functional direction corresponding to the zero frequency mode repre- 
sents an infinitesimal translation of the kink that costs no energy. Since this mode has no restoring 
force, any arbitrarily large amplitude fluctuation along this direction is energetically allowed and 
therefore must be treated non-perturbatively. The variable xq, i.e. the center of mass of the kink 
is elevated to the status of a quantum mechanical variable, and the fluctuations are orthogonal to 
the zero mode. This treatment is the basis of the collective coordinate method[16, [l^, 20]-[^]. 



In collective coordinates quantization instead of the expansion (Q) with (Jsj) we expand <f)(x,t) 



as 



(x,t) = 4> s {x - X (t)) + ^ Qn{t)U n (x - X (t)) 



(10) 
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This amounts to a change of basis in functional space, from the "cartesian" coordinates {q n } 



to "curvilinear" coordinates {xq, Qn^olfl6|i Q| p4 27 1 . 



The next step is to express the Hamiltonian in terms of the new variables xo(t) and Q n (t)- 



For this we follow references [16, [T^, 22, [27|] an d which we summarize below for the cases under 
consideration. 



2.1 Kinetic and potential energies 



In the Schroedinger representation the kinetic energy can be expressed as a functional derivative 
as 



where the functional derivative is written in the new coordinates using the chain-rule 

5 5xq 5 ^ y-, SQ m 5 



x) 5<j){x) 5x ^ q 5(f)(x) SQ r 



(12) 



Taking the functional variation of the field (f>, eq.(10), we obtain 



■r 



S6(x) v-^ 5<f>(x) ^ 

- ^-8xq + }2 -^SQm 



5xo 



m^O 



SQ r 



t> s {X - Xq) 



dx 







dU m (x - Xq) 
Bin 



5x + ^2U n (x - x )5Q n . (13) 



Projecting both sides of the above equation on Uq{x — Xq) and then U*(x — xq) with p / 0, 
using eqn.(p[) and the orthonormalization condition eqn.(p|), we obtain: 



5xp 
5cf){x) 

6<(>(x) 



1 



WJx - z ) 



l + (l/VM)E^o QmS, 

1 En^O GpnQ 



Uq{x - Xq) 



1 + (1/VM) E m / QmS„ 



Uq{x - Xq), 



(14) 
(15) 



where the matrix elements G pm are defined as 



Gpm = J dxU*(x - Xq) 
S m = Gom = / dxU Q (x - Xq 



(X - Xq) 



dx 

dU m {x - Xq) 

dx 



(16) 
(17) 
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At this stage it is straightforward to follow the procedure detailed in|T^, 22, ^7]] to find the final 
form of the kinetic term in the Hamiltonian in the Schroedinger representation of the coordinates 



2\D 5x 5x y/D Sxq 



E 

p,m^0 - 



Gpm Qm 



+ 



Gp m Qm 



y/D 5Q P SQ p x/D 



+ 



E 



D ^ ,„ 5Q P 



-p,q 



VD + ^k Gqn Q n 



D 



5Q„ 



(18) 



where y/D is the Jacobian associated with the change of coordinates [16, 19, 27] and given by 

1 



sf~D = Vm 



1 + 



m^O 



(19) 



The total potential energy, including the elastic term, V[<j>] (see eqn. ([!])), is given by 



V[<f>] = / dx 



1 

2 \ dx 



+ U{4>) 



(20) 



Using the expansion given by eqn.(]i~0|) we find that it can be written in terms of the new 
coordinates as 

V[4>] =M+lj2 QmQ-m^m + 0(Q 3 ) + • ■ • . (21) 

By translational invariance the potential energy does not depend on the collective coordinate. 
Identifying the canonical momenta conjugate to xq, Q n as 



dxo 



for k + . 



(22) 



and using the commutation relation of \/D and 1/y/D with Q n , ir n and P given by 



-iS„ and 



1 



7T n 



D ' 



(23) 



we find the final form of the Hamiltonian: 
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H 




GpmQm , GpmQ 

— 7T-p + TT-p— 5 - 



D 



D 



5-p,q\^D + — P __ — Gq n Q n 



D 



7T_ 



+ X) u p QpQ-p + 

p^O 

-0(Q 3 ) + ---, 



(24) 



where Q„ are now operators. The coordinates Qk associated with the scattering states describe 



the meson degrees of freedom with frequency uj (k) 



u 2 Q = d'U^/d 2 ^. Since 



the Hamiltonian does not depend on xq its canonical momentum P is conserved, it is identified 
with the total momentum of the kink-meson system[16, [n], p2[ . The kink velocity, however, is not 
proportional to P and depends on the momentum of the meson field. 

Since our goal is to study the dynamics of the kink by obtaining the equation of motion for the 
expectation value of the kink collective coordinate, we introduce an external source term linearly 
coupled to xq. This source term has a dual purpose, one is to allow to obtain the correlation 
function of the collective coordinate as functional derivatives with respect to this source, the other 
is to use this source as a Lagrange multiplier to turn the evolution equation into an initial value 
problem. This second use will become clear later when we study the solutions to the equations of 
motion. Therefore we add the term j(t)xo to the Hamiltonian. 



3 A DOMAIN WALL IN THE MESON HEAT BATH 



Our goal is to study the dynamics of a domain wall in interaction with the mesons. This is 
achieved by obtaining the real-time equations of motion of the collective coordinate x$ by treating 



the mesons as a "bath" and obtaining an influence functional[29j-[34] by "tracing out" the meson 
degrees of freedom and the excited states of the kink. We assume that the total density matrix for 
the kink-meson system decouples at the initial time ti, i.e. 



p(ti) = p s {ti) (8) p R (ti 



(25) 



where p s (ti) is the density matrix of the system which is taken to be that of a free particle associated 
with the collective coordinate of the kink, i.e. p s {ti) = \xq >< xq\ and PR(ti) is the density matrix 
of the meson bath and describes mesons in thermal equilibrium at a temperature T. 

Since the kinks can never be separated from the meson fluctuations, this factorization must be 
understood to hold in the limit in which the initial time ij — > — oo with an adiabatic switching of 
the kink-meson interaction. 

The time evolution is completely contained in the time dependent density matrix 

p(t) = U(t,t i )p(t l )U- 1 (t,t i ) (26) 
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with U(t, tj) the time evolution operator. Real time non-equilibrium expectation values and corre- 
lation functions can be obtained via functional derivatives with respect to sources of the generating 
functional ||]-|2| 



Z[j + ,r\ =TrU(oo,-oo;j + )p l U-\oo,-oo;j-)/Trp(t l ), (27) 

where are sources coupled to the meson field and the collective coordinate. This generating 
functional is readily obtained using the Schwinger-Keldysh method which involves a path integral 



in a complex contour in time|35]-[43|: a branch corresponding to the time evolution forward, a 
backward branch corresponding to the inverse time evolution operator and a branch along the 
imaginary time axis from ti to t{ — i/3 to represent the initial thermal density matrix. We will 
obtain the equation of motion for the kink collective coordinate in an expansion of the "adiabatic" 
parameter u;o/Mwhich is also the weak coupling limit of the scalar field theories under consideration 
|16|| . As it will be shown explicitly below in the particular cases studied, the matrix elements given 
by eqns. (|l6| , |r^ ) will provide the necessary powers of the meson mass ujq. The lowest order in 
ijOq/M is formally obtained by keeping only the 1/M terms in the Hamiltonian and neglecting the 
non-linear 0(Q 5 ) terms. Under these approximations, 1/D ~ 1/M and the Hamiltonian has the 
following form 



H = M + 



1 

2M 



where we define 



ran — K -'—mri' 



7T_ m + LJ^QmQ-m + j(t)x , (28) 



(29) 



At this point it proves convenient to write the coordinates and momenta of the mesons in terms 
of creation and annihilation operators obeying the standard Bose commutation relations, 

1 r + i . [uJk 

~1Y 



Qk 



afc + a 



ak- a 



(30) 



The Hamiltonian can be expressed in terms of a and a) as 

H= 2m( P + F[at ' °0 + ^ Uk (°* ak + 1/2 ) + + M ' 



where 



F[a\a] = [ T ik ( a p a k ~ a -p a -k) + Tpk 1 ( a - P a k 

p,fc^0 



(31) 
(32) 



IS) 

We have made use of the symmetries of the operators and defined the symmetric T^ k and 

D kp : :/;;,/'=- J" +\/— A,/,- (33) 



antisymmetric matrices that provide the interaction vertices as 



T (S) = J_ 




j ^k _|_ Up 
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To use the path integral formulation we need the Lagrangian, which to the order that we are 
working {0{ojq/M)) and properly accounting for normal ordering, is given by 



C[xq, a, a" 1 "] = ^-x 2 - x F[(J, a] - ^ u: k (a\a k + 1/2) - j(t)x - M. (34) 

k^O 



The interaction of the collective coordinate and the mesons is now clear. Only time derivatives 
of the collective coordinate couple, a consequence of the Goldstone character of the collective 
coordinate. There are two processes described by the interaction: i) creation and destruction of 
two mesons and ii) scattering of mesons. Whereas the first type can contribute with the mesons 
in their ground state, the second can only contribute if the meson states are occupied. The two 
processes are depicted in fig.(^). As it will become clear below, the second type of processes will 
lead to Landau damping. 

Since we have preferred to work in terms of the creation and annihilation operators it is conve- 
nient to write the path integral for the non-equilibrium generating functional in the coherent state 
representation^, Q. 



Following the steps outlined in reference[13, 44] we find the generating functional of non- 



equilibrium Green's functions in the coherent state representation to be given by 



Z[j + ,j~] = J Vx + J Vx~ J £> 2 7 + J V 2 1 ~exp\i J dt (C[± + , 



7* + ,7 + ,J + ] " C[x-,>y-, 7 *-,r])j 

(35) 



with the Lagrangian density defined on each branch given by 



£[±±,7^,7*^] 



x ± F[i 



dt 



(36) 



and with proper boundary conditions on the fields that reflect the factorized initial condition with 
the mesons in thermal equilibrium. The signs ± in the above expressions correspond to the fields 
and sources on the forward (+) and backward (— ) branches. The contribution from the branch 
along the imaginary time is cancelled by the normalization factor. This is the non-equilibrium 
generalization of the coherent state path integrals (see details in appendix |A|). Non-equilibrium 
Green's functions are now obtained as functional derivatives with respect to the sources . There 
are 4 types of free meson propagators [35]- [42]: 

<4 + (t)a+(0 > = Sk^e-^'^W -t) + n k ] 
< a+(r)at+(r') > = 5 kj> e^-^[9(t - t') + n k ] 
<4 (±) (*)a{(±)(0> = 
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<a { t\t)a^(t')> = (37) 

<4 + (r)a p (r')> = S^e'^^l + n k ] 

<o+(r)o{-(7 y )> = 4, P e^(^')n fc , 

where n k is Bose Einstein distribution for mesons of quantum number k and < • • • > refer to 

averages in the initial density matrix. The ++ ( ) propagators correspond to the time-ordered 

(anti-time-ordered), whereas the ±=p are the Wightman functions. 

An important point to notice is that 

< F[a\a] > = (38) 
in the non-interacting case, since it is proportional to J2k D k -k = 0. 



3.1 The equation of motion for the collective coordinate 

The equation of motion of the expectation value of the collective coordinate for the kink < xq >= 
q(t), can be derived by expanding x (t) = q(t) + £ (t) and requiring <£(£)> = to all orders 
in perturbation theory p6fl . Imposing the condition < £ + (t') >= 0, treating the interaction term 
up to second order in perturbation theory and using eqn.(|38|), we obtain the following linearized 
equation of motion 



dt' < £ + (t)£+(0 > \{Mq(t') + f dt"T m (t' - t")q(t")\ + < tm + (t') > j(t'] 



o.. 

(39) 



where the retarded kernel is given by 

-iT m {t-t')6{t-t') = <F[a^ + (t),a + {t)}F[a^ + (t'),a + (t')}> 

- <F[a^ + (t),a + {t)}F[a^-(t').a-{t')}> (40) 

Since we restrict ourselves to non-relativistic kinks we consider g < 1, The non-equilibrium 
Feynman diagrams that contribute to one loop order (second order) are shown in fig.(^|). 



Alternatively this equation of motion may be obtained by computing the influence functional [29] - 
[34] in second order perturbation theory. The resulting influence functional is quadratic in the 
collective coordinate, performing the shift x^(t) = q(t) + £ (£) the above equation of motion is 
obtained by requesting that the linear terms in £ vanish, (there are two linear terms, both give 
the same equation of motion). 

The kernel T m (t — t') is found by using eqn. © and eqn.([3|) and it is given by 
T m {t-t') = -4 {T$T[ll k (l+2n p ) S m[(u p + u; k )(t-t')] 

- 2T^ ) T[f_ k n p sm[(u Jp -u Jk )(t-t')]}. (41) 
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Performing the integral over t' in eqn.(|3^) by parts, we obtain the final form of the equation of 
motion ^ 

Mq(t) + / dt' £ m (t - t')q{t') = j(t), (42) 



— oo 



where the non-local kernel is given by 

^_ dT m (t-t') _ dT m (t-t') 

^m(t-t)- - - — . (43) 

Using eqn.([il]) we find the final expression for the kernel S m : 

£ m (t-t') = _ 4 J2 {T p ( f ) ^^ ) _ fe (l + 2n p )(w p + w fc )cos[(a; p + w fc )(^-^ , )] 



p,k^0 



(A) rp[A) 



-2 T£' T 



f-k n P ( w p - cos - Uk){t ~ t')} } . (44) 



We will see in the next sections that the two kernels S m ; T m have very special significance: whereas 
S m is identified with the real-time retarded self-energy of the collective coordinate, T m will provide 
the coefficient of dynamical friction in the Markovian approximation. 

It is more convenient to express the equation of motion of the kink in terms of the velocity 

MV(t) + f dt' £ m (t - V(t') = j(t) (45) 

with E m given by eqn.(E3). 



The relation (|44|), ensures to this order in the perturbative expansion, that with an adiabatic 
switching on convergence factor introduced to regularize the lower limit of the integral and to 
provide an initial factorization of the density matrix as t{ — > — oo the total integral of the retarded 
self-energy kernel vanishes, i.e, 

rt 

T, m (t - t')alt' = 0. (46) 



This result is consistent with that found in Refs. [14, 15 1. 



Therefore for j = 0, any constant velocity of the kink is a solution of the equation of motion 
|). This result is physically clear: when j = 0, the total Hamiltonian commutes with P, the 
canonical momentum conjugate to xq because of translational invariance, i.e. the total momentum 
of the system is conserved. One can then go to a frame in which P = and since the meson bath 
is in equilibrium this must result in that the domain wall must have a constant velocity, therefore 
for j = there must be a constant velocity solution to the equations of motion of the collective 
coordinate resulting in 
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3.2 General properties of the solution 



Since in the absence of an external driving term we have found that the domain wall moves with 
constant velocity, we can now use the external source term to cast the evolution as an initial value 
problem. For this consider the situation in which at time t = a force is applied, therefore changing 
the velocity of the domain wall. Assuming that for t < the kink traveled with a constant velocity 
Vq, after switching on the external force the domain wall will accelerate, but it will also transfer 
energy and excite the meson degrees of freedom and this will lead to dissipative processes. Therefore 
writing V(t) = v o + v(t) with j(t < 0) = ; j(t > 0) ^ and using the property (p3|) the equation 
of motion for the velocity change becomes 

Mv(t) + [ tit' E m (t - v(t') = j(t). (47) 
Jo 

The solution of this equation is found by Laplace transform, in terms of v(s) ; S m (s) ; J(s), 
the Laplace transforms of the velocity, self-energy kernel and current respectively, in terms of the 
Laplace variable s. We find that the solution is given by 

and consistently with the linearized equation of motion and the non-relativistic dynamics v(t); t>o 
1. The quantity 

G{s) = i iV n (49) 

is the Laplace transform of the propagator of the velocity of the collective coordinate. The real 
time evolution is found by the inverse Laplace transform 

v(t) = — [ e st v{s)ds (50) 
2m Jc 

where C refers to the Bromwich contour running along the imaginary axis to the right of all the 
singularities of v(s) in the complex s plane. Therefore we need to understand the analytic structure 
of G(s) in eqn. ( f48|) to obtain the real time dynamics. The Laplace transform of the self-energy 
kernel is conveniently written as a dispersion relation in the form 



S m (s) = st m (s) (51) 
f m (s) = [ 4r^2 d Po (52) 
p{Po) = -4 {^kT-J-ki 1 +2n p )S(p -u p -u k ) 

p,k^0 

- T { ^_ k {n p - n k )5{p - uj p + u; k )} , (53) 
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where T m (s) is the Laplace transform of the kernel T m given above. 

This dispersive form for the Laplace transform of the kernel reveals that r m (s) has a disconti- 
nuity in the complex s-plane along the imaginary axis, since 

f mI (s = iu + 0±) = [p(M ) _ ^(-1^1)] (54) 

The imaginary part changes sign with w as a result of the retarded nature of the kernel. There- 
fore the propagator G(s) has cuts along the imaginary axis in the complex s-plane. The two different 
contributions to the spectral density (53) yield to two different cut structures. For uj > 0, the first 



term, proportional to 5(u> — ui p — ujk) gives a two-meson cut beginning at 2ujq corresponding to the 
process of spontaneous and induced two-meson creation and annihilation. The second contribution 
corresponding to 8{uj — uj p + uJk) gives a cut, which we identify as the Landau damping cut[^7], 46], 



pinching the origin and originates in the process of scattering of mesons present in the medium off 
the domain wall. As it will be seen in detail for the examples in the next section the discontinuity 
vanishes linearly as uj — > allowing an analytic continuation into the second Riemann sheet and 
to isolate the pole. This linear vanishing of the self energy is consistent with the case studied 



byjl4|, 15]. This process is present only for finite temperature as there must be mesons present 
for this process to exist. This contribution is identified as Landau damping from the in medium 
mesons and will be seen to provide the leading contribution to the long time relaxation. 

The presence of a static friction coefficient will be revealed by a pole in G(s) with a negative 
real part, since this will translate into an exponential relaxation of the velocity. 

In the absence of interactions G(s) has a simple pole at s = 0. Since we obtained the expres- 
sion for the kernels in perturbation theory the position of a pole must be found in a consistent 
perturbative expansion by writing s p = {1/M)s\ +■■■■, we find 

s p = --^£ m (s = 0) = 0. (55) 

Therefore the coefficient of static friction vanishes. This is a consequence of the vanishing of the 
integral ([46|). Therefore up to this order in perturbation theory the position of the pole in the 
s- variable remains at s = resulting in that the static friction coefficient vanishes. 

In summary, the analytic structure of G(s) in the complex s-plane corresponds to a pole at 
s = with residue 

Z s = X (56) 

i + i?r m (o) 

and cuts along the imaginary axis beginning at ±2iuj Q ; ±ie with e —* to clarify that the beginning 
of this cut pinches the pole at the origin but the continuum contribution to the spectral density 
(discontinuity) vanishes at the position of the pole at s = 0. 

The residue Z s has a very clear interpretation, it is the "wave function renormalization" and 
its effect can be understood in two alternative manners. 
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Consider the case in which J = in eqn.(|48|). Performing the inverse Laplace transform and 
invoking the Riemann-Lebesgue lemma, the long time behavior will be completely dominated by 
the pole at s = 0. Therefore, if the velocity of the kink has been changed at t = by some external 
source, this disturbance will relax in time to an asymptotic value given by 

Voo = Z s v . (57) 



Alternatively, consider the case of vq = but with an external source term switched on at 
t = and constant in time thereafter. Again the inverse Laplace transform at long time will be 
dominated by the pole, and we find that the kink moves with constant acceleration given by 



J 



M, 



with M, 



eff 



M 



(58) 



Thus the wave function renormalization can also be understood as a renormalization of the kink 
mass. The ratio of the asymptotic acceleration to the initial acceleration is given by Z s . As the 
kink moves, the interaction with the meson bath "dress" it changing its effective mass, which will 
be seen in specific models to be larger than the bare mass. 

Furthermore we can now derive the following important sum rule. Consider the case j = 0. 
Isolating the contribution from the pole and the cuts (by replacing 5{p —uj p +ujk) — ► 5(p — w p +u;fc— e) 
and taking the limit e — > at the end of the computation) separating the pole and continuum 
contributions, and deforming the contour of integration for the inverse Laplace transform as shown 
in fig. (|^) we find the time evolution (for j = 0) to be given by 



v(t) 



Z s + 



duj 



t Im (u) cos(wt) 



Vq ttM 
Evaluating at t = we obtain the sum rule 



1 + 



M 



+ 



Tim (a;) 
M 



(59) 



Z s + 



ttM 



du 



-T Im ((j) 



Y _|_ ^Rm {w) 



M 



+ 



ri m (tj) 

M 



(60) 



Since the spectral density p{oj) is positive (semi) definite as it will be explicitly shown below for 
specific models, the sum rule above determines that 



Z R < 1 



< 1 



vo 



(61) 



Although a sum rule similar to eqn.(|60|) is obtained in quantum field theory from the canonical 
commutation relations, its validity for the collective coordinate associated with the domain wall is 
far from obvious since the kinematic and canonical momentum conjugate to the collective coordinate 
are different. 
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The continuum contribution in eqn.(59) is dominated at long times by the small u region. 
Therefore for T / the asymptotic long time relaxation of the velocity is completely determined 
by the Landau damping cut which has support at small u>, whereas the two-meson cut vanishes 
below the threshold at 2lvo. 

A further understanding of the dynamics will necessarily require knowledge of the matrix ele- 
ments to establish the details of the kernels. This will be studied in particular models in the next 
section. 



3.3 Semiclassical Langevin equation 

The classical Langevin equation is an adequate phenomenological description of Brownian motion 
obtained by considering the dynamics of one (or few) degrees of freedom that interact with a bath 
in equilibrium. 

It contains a term proportional to the velocity of the particle which incorporates friction and 
dissipation and a stochastic term which reflects the random interaction of the heat bath with 
the particle. These two terms are related by the classical fluctuation-dissipation relation which is 
derived in appendix [C] (see eqn. ( 130 ). 



At the quantum mechanical level it is also possible to obtain a "reduced" or coarse grained 
description of the dynamics of one (or few) degrees of freedom in interaction with a bath. The 
coarse graining procedure has a very precise meaning: the full time dependent density matrix is 
traced over the bath degrees of freedom yielding an effective or "reduced" density matrix for the 
degrees of freedom whose dynamics is studied. 

Such a description of non-equilibrium dynamics of a quantum mechanical particle coupled to a 
dissipative environment by a Langevin equation was presented by Caldeira and Leggett p0| and by 



Schmid [31]. Their technique is based on the influence- functional method of Feynman and Vernon 



[29] that naturally leads to a semiclassical Langevin equation. 



In this section we follow the procedure of f|3C|l -|34| generalized to our case to derive the Langevin 



equation for kinks in a heat bath to lowest order in the adiabatic (weak) coupling. 

The main step is to perform the path integrals over the meson degrees of freedom, thus obtaining 
an effective functional for the collective coordinate of the kink. Unlike the most usually studied 



cases of a particle linearly coupled to an harmonic reservoir] 30 1-| 34] we have here a bilinear coupling 
to the mesons. Therefore the influence functional cannot be obtained exactly, but it can be obtained 
in a consistent perturbative expansion. For this we treat the interaction term Ci[x^, ^ , 7* 1 * 1 ] in 
perturbation theory up to second order in the vertex proportional to x^ (which is equivalent to 
lowest order in the adiabatic coupling m/M). Integrating over the meson variables and using 
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< F[aJ, a] >= 0, we obtain 

2[j + ,r = o] 

where 



Vx 



+ Vx ~ e ir.<^[* + ]-£o[* i) ^-[i+.i-], 



±\2 • ± 



£0^] = -M^Y-j 



(62) 
(63) 



and ^"[x+ji; ] is the influence functional [ 29 ]- ■ To lowest adiabatic order we find 



J : [x + , x ] = exp I — — J dtdt' 



x + {t)G ++ (t,t')x+(t') + ±-(t)G—(t,t')x-(t') 



+ x + (t)G + -(t,t')x~(t') + ±-{t)G- + (t,t')x + (t ') 

in terms of the real-time meson correlation functions (see appendix (^)) 

G ++ (t,t') = <F[J + {t),a + (t)}F[J+(t'),a + {t')}> 

G—(t,t') = < F[(J-(t),ar(t)] F[a}-(t'),cr(t')]> 

G + -(t,t') = - <F[a^ + (t),a + {t)}F[a^-(t'),a-(t')}> 

G- + {t,t') = - < F[<f-(t),ar(t)] F[a^ + (t'),a + (t')} > . 



(64) 



(65) 



At this stage it is convenient to introduce the center of mass and relative coordinates, x and R 
respectively, which are defined as 



R(t) = x + (t) - x~(t). 



(66) 



These are recognized as the coordinates used in the Wigner transform of the density matrix 1 3C ]- [ 34 ] 
in terms of which the partition function becomes 



(67) 



Z[0] = JVxVR e iS[x > R] 
with the non-equilibrium effective action given by 

S[x,R] = J dtR{t) -Mx(t) - ^ dt' (^(t - t') ±(t') - K(t-t')R{t') 
in terms of the kernels K\(t — t') and K(t — t') which are given by (see appendix (^) ) 

K^t-t') = 8*0(t-O \r^ ] T[ll k (l +n p + n k ){u p + u k ) cos [(u p + u k )(t-i/) 



(68) 



p,fc^0 

<A) rp(A) 



' T pk T -p-k( n P - n k) i^P ~ ^k) cos [( W P ~ u k)(t ~ t')] 

-2iE m (t-0 



(69) 
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and 



Kit - i!) 



-2 \ T pk^ T -p-fc( 1 +n p + n k + n p n k )(uj p + u k ) 2 cos [{u p + u k )(t - t') 
+ 2 T^_ h n k {l + n p ) (w p - w fc ) 2 cos [{u p - oj k )(t - t')] X. 



(70) 



At this stage it proves convenient to introduce the identity 

e -i/ dtdt'R(t) K(t-t')R(t') _ c ^ f V £ e -ydtdt'£(t)K-Ht-*)W)+ifdt£(t)R(t) ^ 



with C(t) being an inessential normalization factor, to cast the non-equilibrium effective action 
of the collective coordinate in terms of a stochastic noise variable with a definite probability 



distribution! 31 1- 1 34]. 



Mx(t) / d£ Ki(t - t') x(t') + £(t) 



Z[0] = J VxVRV£P[£\ expS^i J dtR(t) 



where the probability distribution of the stochastic noise, P[£], is given by 

P[£] = J Vtexp{-^j dtdt't(t)K-\t-t')ti(t' 



(72) 



(73) 



In this approximation we find that the noise is Gaussian, additive and with correlation function 
given by 

<mm>= K(t-t>). (74) 



respect to E(t)l3q]-|M 



The semiclassical Langevin equation is obtained by extremizing the effective action in eqn. (phwith 
Mxit) + f dt'Z m {t-t')x{t')-j{t)=£{t). (75) 



Two features of the semiclassical Langevin equation deserve comment. The first is that the 
kernel K\{t — t'), as can be seen from eqn.(|69|), is non-Markovian. The second is that the noise 
correlation function K{t — tf) is colored, i.e. it is not a delta function S(t — tf). The relationship 
between the kernels K\(t — t') and K(t — t') established in appendix (|C]) constitutes a generalized 



quantum fluctuation dissipation relation pC|]-[ 34] (see eqn.( |130] ). Finally we recognize that taking 
the average of (|75| ) with the noise probability distribution P[£] yields the equation of motion for 
the expectation value of the collective coordinate (eqn. ([45])). 
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A classical description is expected to emerge when the occupation distribution for the mesons 
can be approximated by their classical counterparts |Q, i.e. when rik ~ T/u>k- 

//the kernels T, m and K admit a Markovian limit then a diffusion coefficient could be extracted 
by computing the long time limit of the correlation function << (x(t) — x(0)) 2 >> /t where 
<< • • • >> stand for average over the noise distribution function. However when the kernels do 
not become Markovian, such a definition is not appropriate. 

This summarizes the general formulation of the description of the dynamics of the collective 
coordinate both at the level of the evolution equation for the expectation value as well as for the 
effective Langevin dynamics in terms of stochastic noise terms arising from the fluctuations in the 
meson bath. We are now in condition to study specific models. 



4 TWO MODELS 



In the previous sections we established the general aspects of the real-time dynamics of kinks in 
the presence of the meson bath, obtaining the equation of motion as well as the Langevin equation 
for the collective coordinate in lowest adiabatic order. Further progress in the understanding of the 
dynamics necessarily involves the details of particular models which determine the matrix elements 
ji(A,s) an( j th ere f ore th e time dependence of the kernels involved. In this section we study these 
details for the Sine-Gordon and </> 4 models. 



4.1 Sine-Gordon 



For the Sine-Gordon model the potential is given by 

A 



COS 



V>4 



and the static kink solution is given by|L7|, [16, [26[] 

4 



4> s (x) = — ^=arctan |e'' 
V A 



the kink mass and the adiabatic ratio are given by 

8m m A 
~T ' M ~ 8 



M 



(76) 



(77) 



(78) 



The normal modes of this theory are the solutions of the following equation, (see eqn.(|6])) 



2m 1 



dx 2 „„„u2 



cosh (mi) 



lp n (x) = Unlp n (x). 



(79) 
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The solutions of the above differential equation are well known p5|| , |r7|| , pq| . There is only one 
bound state with zero eigenvalue, the zero mode, followed by a continuum with wave functions 
given by 



U k (x) 



1 



2lTUJ k 



(—ik + mtanh(mx))e 



ikx 



(80) 



with lo\ = k 2 + m 2 , i.e. wq = m. The scattering states represent the meson excitations around the 

kmkra. 



The matrix elements D pk were already calculated by de Vega|2q], (see eqns.(2£,16|)) and are 
given by 

i(p 2 — k 2 ) 



D kp = ip5(k+p) + 



TT (p+k) 

2 m 



Auj k u: v sinh 

which determine the symmetric and antisymmetric matrix elements 



for p 7^ k 



(81) 



T (S) 

pi 



T (A) 



UJ, 



1/2 



1/2 



1/2' 



1/2' 





o 2 ) 


AiOk^Op sinh 


'it (q+p)~ 

2 m 



— I + — 



UJr, 



CO*, 



Auj k uj p sinh ^ 



tt (g+p) 

2 m 



(82) 



Since in this theory there are no bound states other than the zero mode, F[a^ , a] is given only 
by the first two terms in eqn.( |12 j ). We recognize the "structure factor" 



S(Q) 



dx 



2m 2 



Q 



(83) 



, 2vr cosh 2 [mx] sinh[§^] 

This structure factor will play an important role in understanding the large energy behavior of the 
one-loop contribution. The important point to notice is that the structure factor is dominated by 
momenta Q ~ m, falling off exponentially for \Q\ 3> m. 



case 



Substituting eqn.(82) in eqns.(|4l|) and eqn.(|44"|), we obtain the final form of the kernels in this 
] /' x ,, / :x „ (Q - 2k) 2 



fm (t t 



4 3 



dQ\S{Q)Y 



dk- 



k^k 



|(1 + 2n k )(ujQ_ k - w fc ) 2 sin [(ujQ-k + u k ){t - t')] 



(n Q „ k - n k )(ujQ__ k + uj k ) 2 sin [{u Q - k - UJ h )(t - t')] } 



(84) 



1 C°° C°° (O — 2k) 3 
S m (t-0 = -3/ dQ\S{Q)\ 2 Q dk { ^ H(l + ng_fc + n fc )(w g _ fc -w fc )x 

4 J — 00 J — OO ^Q— k ^ k 

cos [(wQ_ fc + w fc )(t -i')] - (n Q _ fc -n fc )(a;Q_ fc + a; fc )cos [(w Q „ fc -w fc )(t-t')]} • 

(85) 
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The introduction of S(Q) clarifies that Q is the momentum transferred into the meson loop, 
and because S(Q) is peaked at Q = with a width of the order of the meson mass we conclude that 
the momentum transferred into the meson loop is of the order of the meson mass. This observation 
will prove to be very important in the analysis of the high temperature limit in a later section. 

It proves useful to express T^ m (t — £') and E m (£ — £') in terms of dimensionless Quantities to 
display at once the nature of the adiabatic expansion. To achieve this let us make the following 
change of variables 

Q k T 

Q — > — ; k — > — ; t = mt and T = — . (86) 
mm m 

Then T m (t — t') and E m (t — t') can be written as 

T m (t - 1') = m 2 T(r - t') and E m (t - t') = m 3 E(T - r'), (87) 

where 

/oo 
dQdk\S(Q)\ 2 {Ti(Q, k) sin [(w Q . k + w k )(r)] + T 2 (Q, k) sin [(w Q „ k - w k )(r)}} (88) 



— oo 
oo 



S(r) = / dQdk\S(Q)\ 2 {Si(Q, k) cos [(w Q ^ k + w k ){r)) + E 2 (Q, fc) cos [(w Q _ k - w k ){r))} 
with 

r fn u\ 1 (l + 2n fc )(Q-2fc) 2 ( WQ _ fc - w k ) 2 

l(Q ' } " 64 

= J_n fc (Q-2)b)> q fc + Wfc )^ 

1 j 32 w%_ k wl 

Ei(Q,k) = (wQ- k + w k ) ri(Q,fe) 

S2(p,fc) = (i«q-a; - ^fc) r 2 (Q,fe) 

= A: 2 + 1 ; n fc = . (90) 

e r — 1 

Fig.(^J) shows the numerical evaluation of T(r) and S(r) vs. r for different values of T. We 
clearly see that the self-energy kernel E is peaked near r = and localized within a time scale 
t s ~ Similarly, the kernel T varies slowly over a large time scale ~ 5 — 10m _1 . 



4.1.1 Equation of motion: Exact solution vs. Markovian approximation 

In terms of dimensionless quantities the equation of motion (f45|) becomes in this case 



v(t) + - / T dr ' S ( r - t')v(t') = J, (91) 
8 Jo 



20 



where J = j / (mM) and the dot stands for derivative with respect to the dimensionless variable r. 

As shown in fig.(|Cj), the kernel S(r) has "memory" on time scales a few times the inverse of the 
meson frequency. If the velocity of the domain wall varies on time scales larger than the "memory" 
of the kernel a Markovian approximation to the dynamics may be reasonable. The first step in the 
Markovian approximation corresponds to replacing v(t') by v(t) inside the integral in eqn. ( |9l| ) 
and taking it outside the integral. A second stage of approximation would take the upper limit of 
the integral to oo thus integrating the peak of the kernel. However we have shown above that the 
total integral of the kernel vanishes, thus this second stage cannot be invoked. Recognizing that 
Jq S(t — T / )dT / = r(r) the Markovian approximation to (|9l| ) is given by 

v(t) + ~«(T)r(r) = J. (92) 

As advanced in the previous section, we now identify the kernel T(r) as the dynamical friction 
coefficient in the Markovian approximation. The property ( f46[) determines that T(t — > oo) = 0. 

We will now focus on the initial value problem with J = and v(t = 0) = vq- The formal 
solution of the equation of motion in the Markovian approximation is given by 

v(r)=v e-^o r ^ (93) 

Even in the Markovian approximation the relaxation of the velocity at long times is not exponential 
because T(r) — > at long times as can be seen in fig. (|C]). 



4.1.2 Velocity relaxation and wave function renormalization 

In order to display more clearly the dissipative effects, we now study the relaxation of the kink 
velocity. For this consider the initial value problem with j(t > 0) = and initial velocity v{t = 
0) = v . 

As the kink moves through the bath, its velocity decreases because of the interaction with the 
mesons, the asymptotic final velocity is related to the initial velocity through the wave function 
renormalization as explained section |3.2| above. We present the numerical solution of the homoge- 
neous equation for v(t)/vo in figure (|C]), where we also present the homogeneous solution in the 
Markovian approximation described above. We clearly see that the initial velocity relaxes to an 
asymptotic value Voo/vq. However the time dependence cannot be fit with an exponential. We 
can see that even at high temperatures the Markovian approximation grossly fails to describe the 
dynamics. 

According to the analysis of the general solution, the ratio v^/vq should be given by the wave 
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function re-normalization, i.e. 

Z = - = — (94) 

S l + $f(s = 0) V 

Table 1 below compares the ratio foo/^o obtained from the numerical solution to the exact 
evolution equation, with the value of the wave-function renormalization. Clearly the agreement is 
excellent, confirming the analysis of the asymptotic behavior of the solution in real time. 



Table (1): Numerical evaluation of Z s and Voo/vq in Sine-Gordon theory for m/M = 0.1,0.25. 







z s 




m/M = 0.1 


m/M = 0.25 


m/M = 0.1 


m/M = 0.25 


Zero Temp. 


0.999808 


0.999521 


0.999808 


0.999521 


Temp. 1.0 


0.993438 


0.983754 


0.993438 


0.983753 


Temp. 5.0 


0.96055 


0.906885 


0.96055 


0.90687 


Temp. 10.0 


0.923458 


0.828352 


0.923446 


0.828303 



4.1.3 Kernels for the semiclassical Langevin equation 



Knowledge of the matrix elements T^ A > , allow us to obtain the final form of the kernels that 
enter in the semiclassical Langevin equation given by eqns.(|69|) and (|70|), and eqn.(|82[). These 
kernels can be written in terms of the dimensionless quantities given by eqn.(|86[). Since K\(t — t') = 
) we focus on K(t — t'). In term of dimensionless quantities, K (t) = m 4 /C(r) where 



/C(r) = / dQdkQ 2 \S(Q)\ 2 {d{Q,k) cos [{w Q _ k + w k )(r)) + C 2 (Q,k) cos [(w Q „ k - w k ){r)]} . 

(95) 



with 



= 1 (l + "Q-* + "* + "Q - k n k )(Q-2 k y 

4 4 w Q-k w k 

= l ^UylW^ (97) 

4 W "-k W k 



The contribution from C\ is recognized to arise from the process of emission and annihilation 
(spontaneous and induced) of two mesons, whereas that from C2 arises from the scattering off in 
medium mesons and has its origin in the Landau damping diagram shown in figure (Q). 

Fig.(^|) shows JC(t) for different temperatures T. Notice that at large temperatures the kernel 
becomes strongly peaked at r = and one would be tempted to conclude that the classical limit 
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corresponds to a delta function. However the coefficients (9^,9^) are such that the total integral 



in r (leading to delta functions of sums and differences of frequencies) vanishes. We then conclude 
that even in the high temperature limit the noise-noise correlation function is not a delta function, 
i.e. the noise is "colored" , the classical fluctuation dissipation relation in terms of a delta function 
noise correlation does not emerge and a diffusion coefficient cannot be appropriately defined. We 
postpone until a later section a discussion of the high temperature limit and the classical regime. 



4.2 4 Theory 



In this model the potential is given by 

The static kink solution is given by 

(f> s (x - x ) 

and the kink mass is given by 



m 
2A 



(l-A0 2 ) : 



M 



tanh [m(x — xq)] 



Am 
3A 



(98) 

(99) 
(100) 



and the normal modes are the solutions to the equation, (see eqn.(y)) 

72 Q m 2 



-3 + 4m 2 -j- 

ax z cosh (mx) 



1p n (x) = UJ n 1p n (x). 



(101) 



The solution of the above differential equation is well known [45], [17]. It has two bound states 
followed by a continuum. The normalized eigenvectors are given by 



U (x) 
U b (x) 
U k {x) 



'3m , 2 



sech \mx] oc — - — 
2 L J dx 



with loq = 



'3 m 



sech[mx] tanh[mx] with u> 2 = 3 m 2 



m 2 e lkx f l2r i „ . k , r . 1 k 2 

— — < 3tanh [mx\ — Si — tanhlmxj — l =■ 

y / 2ir(k 2 + m 2 )ujk I m m l 



(102) 



with uj 2 = k 2 +Am 2 . The scattering states are identified with meson modes and the meson frequency 
is identified with uj q = 2m. 

The bound state with zero frequency is the "zero mode" , whereas the bound state with ui 2 = 3m 2 
corresponds to an amplitude distortion |l8|, [ill or excited state of the kink. 



The matrix elements -D p /% are given by, (see eqns.(^|T 
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Dbk 
D p k 



sech 



2m 



Vk 2 + m 2 (k 2 + 3m 2 ) (from the bound state) 



m 2 ojk 



ik5(p + k) + 



3z7r(fc 2 — p 2 ) (p 2 + k 2 + Am 2 



Am 4 N p N k sinh 



7T (p+k) 

2 m 



for p ^ k 



where Nk is defined as 



N k 



'2irwl(k 2 + m 2 



(103) 



(104) 



We notice that the coupling to the continuum through the bound state given by the matrix 
element D b k is of the same order as the coupling to the continuum-continuum (matrix elements 
Dpk). This will have interesting consequences for the dissipational dynamics. The symmetric and 
antisymmetric matrix elements for the continuum states are given by 



T (S) 
VI 



T (A) 
VI 



3_ 

32 

3 

32 



1/2 



\ 1/2 

v\ 



+ 



<jJa 



OJr, 



U>r, 



1/2 



1/2' 



(q 2 -p 2 ){p 2 + q 2 + 4m 2 



\J q 2 + m 2 \Jp 2 + m 2 Uq(jj p sinh 
(q 2 — p 2 )(p 2 + q 2 + 4m 2 



7T (£+P) 

2 m 



\J q 2 + m 2 y/p 2 + m 2 uj q oj p sinh 



7T (g+p) 

2 m 



(105) 



whereas those involving the bound state are obtained by replacing the matrix elements Dbk for the 
Dpk- 

Since in this model there is one bound state other than the zero mode, the interaction vertex 
F\a\ a] is given by eqn.( |123[ ) in the appendix. The contributions from bound-state-continuum vir- 
tual transitions do not mix with the continuum-continuum to this order in the adiabatic expansion. 
As a consequence of this simplification the dimensionless kernels (in terms of the dimensionless 
variables introduced in (j36|)) become 



T(r) = r dpi rj(p) sin [(w p + w b )(r)] + T b 2 (p) sin [(w p - w b )(r)] 

J -co I 

/oo 
dQ|S(Q)| 2 dfc{ri(Q,A:) sm[(w Q -k + w k )(T)} + T 2 (Q,k) sin [{w Q ^ k - w k ){T)]} 
'CO 

S(r) = [ dp\ T,\(p) cos[(wp + w b )(r)] + S|(p) cos [(w p - w b )(r)] 

+ I dQ\S(Q)\ 2 dk{T 11 {Q,k) cos [(u>Q_ fe + w k ){r)] + S 2 (Q, fc) cos [(wQ_ fc - u?fc)(r)]} 



(106) 
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with 



Ti(Q,k) 

r 2 (Q,k) 
it(p) 

Si(Q,fc) 

s 6 i(p) 



32 (1 + + nfcKQ - 2k) 2 (w Q _ k - w k ) 2 ((Q - fc) 2 + fc 2 + 4) 

4 4 ^_ fe <4((Q-fc) 2 + l)(fc 2 + l) 

3^ (n k - n Q ^ k )(Q - 2k) 2 {w Q . k + w k ) 2 ((Q - A:) 2 + fc 2 + 4) 2 
4 4 tW 3ffC?-fc) 2 + l)fA; 2 + l) 



iry/3 (p 4 + 4p 2 + 3) 2 
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[w p - w b )(l + n b + ra p) „ cch 2 



7rp 



7rp 



ttV3 (p 4 + 4p 2 + 3) 2 (w p + w b )(n b - n p ) 2 
128 — itffc) 

(^Q-fc + w^jfc) ri(Q,fc) ; E 2 (Q,k) = (wQ- k - w k ) T 2 {Q,k) 
{w p + w b ) V\(p) ; Sa(p) = Op - w 6 ) r|(p) 



(107) 



where S(r) and T(r) are defined as in eqn.(p^). The functions X(t) and T(r) where evaluated 
numerically at different temperatures T, the results are displayed in fig.(|C|). The behavior of these 
functions differ from those in the Sine-Gordon theory because of the presence of the bound state 
which is interpreted as an excited state of the kink. As the kink moves in the dissipative medium, 
energy is transferred between the kink and the bound state resulting in the Rabi-like oscillations 
displayed in the figure. We notice that the contribution of the bound state is of the same order of 
magnitude as that of the continuum. 



4.2.1 Equation of motion: Exact solution vs. Markovian approximation 

The solution to the equation of motion and the comparison to the Markovian approximation pro- 
ceeds just as in the the case of the Sine-Gordon model. The equation of motion is again solved as 
an initial value problem. The exact and Markovian solutions are displayed in figure (|C]). 

The new feature of the solution are the oscillations that result from virtual transitions to the 
bound state. We interpret these in the following manner: as the kink moves it excites the bound 
state that corresponds to a kink distortion j^], this excitation in turn reacts-back in the dynamics 
of the collective coordinate in a retarded manner. 

While the exact solution in this model is qualitatively similar to that of the Sine-Gordon model, 
we see however, that quantitatively they are different: there is stronger dynamical dissipation in 
the (j> model as compared to the Sine Gordon case, due to the strong coupling to the bound 
state-continuum intermediate states. 
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4.2.2 Velocity relaxation and wave function renormalization 



In this model the Laplace transform of the functions T(t) and S(r) are given by 
E00 = / dQ\S{Q)W 5 El(g ' fc)S ^ + MQ,k)S 



Us) 




S 2 + (iVQ-k + Wk) 2 s 2 + {wQ- k - w k ) 
X\(p)s + X b 2 (p)s \ 



2 



s 2 + (w p + w h ) 2 s 2 + (w p - w b Y z 



oo 



r( s ) = / d Q \s(Q)\ 2 dk \ T l [Q \ k) {Wp + ^ + r2(g ' fe) K " Wk) 



s 2 + (^Q-fc + w k ) 2 s 2 + (w Q _ k - w k ) 2 



+ / dH r;(rtK + 11 ,,) r|( P )K-^) l 

' ' S 2 + (w p + W fe ) 2 S 2 + (w p - W b ) 2 J 



— oo 



With the quantities S 6 ; T b given above. The homogeneous equations of motion given by 
)l]) (exact) and its Markovian approximation ( |92|) both with j = are solved with the kernels 
£(t) ; r(r) given above for v(t)/vQ. The asymptotic behavior of the exact solution will be compared 
with the prediction Voq/vq = Z s , with the wave function renormalization Z s given by eqn.(|94"|) but 
with the r(s = 0) appropriate to the (ft 4 model. 

Figure shows the numerical solutions of eqn.(91) and eqn.(|92|) with j = for (v(t)/vo) 
for temperatures T = 0, 1.0, 5.0 and 10.0. Again the Rabi-like oscillations associated with the 
excitation of the bound state is apparent in the solutions. We have checked numerically that 
asymptotically the velocity tends to a constant value Voo but not exponentially. Table 2 shows the 
values of Voo/vq and Z s for these temperatures for m/M = 0.1 and 0.25 where foo/^o was evaluated 
at r = 200 for the exact solution. Within our numerical errors, we can see that eqn.(|94|) is fulfilled. 



Table (2): Numerical evaluation of Z s and Voq/vq in <ft 4 theory for m/M = 0.1,0.25. 





Voo/Vo 


Z s 




m/M = 0.1 


m/M = 0.25 


m/M = 0.1 


m/M = 0.25 


Zero Temp. 


0.999225 


0.998064 


0.999176 


0.997943 


Temp. 1.0 


0.961561 


0.911004 


0.96376 


0.914072 


Temp. 5.0 


0.784934 


0.593058 


0.787734 


0.597494 


Temp. 10.0 


0.642698 


0.417231 


0.646577 


0.422562 
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4.2.3 Kernels for the semiclassical Langevin equation 



From the definition of the kernels K\(t — t') and K(t — t'), eqns.(39) and (|70|), and eqn.( |105l) , these 
kernels can be written in terms of the dimensionless quantities given by eqn.(p6|) as 



/C^r-r') = -2iS(r-r'), 
where E(r — r') is given by eqn. (|l06| ) and K(t) = m 4 /C(r) with 

/C(r) = 



(109) 



dp { Cl{p) cos l(w p + w b )(r)] + C b 2 {p) cos [(w p - w&)(t)]} 



+ 



dQ|5(Q)| : 



d/c{Ci(Q, fe) cos [(wq_a,. + w fc )(T)] + C 2 (Q,k) cos [(wQ-k - w k )(T)]} 

(110) 

with the dimensionless matrix elements 

18 Q 2 (l + n Q _ k + n k + n Q _ k n k ){Q - 2k) 4 ((Q - k) 2 + k 2 + 4) 2 



Ci(g,fc) 
c b 2 (p) 



4 5 u;3_ fcU; 3 (( g_ fc) 2 + 1)(A; 2 + 1) 

_9_ n fc Q 2 (l + n Q . k ){Q - 2fc) 4 ((Q - k) 2 + k 2 + 4) 2 
4 4 ^_ fe ^((Q-fc) 2 + l)(fc 2 + l) 

ttv/3 (p 4 + 4p 2 + 3) 2 (p 2 + 1)(1 + n b + n p + n b n p ) 
4 4 



w^cosh 2 [^] 



27rV3 (p 4 + Ap 2 + 3) 2 (p 2 + l)ra p (l + 



nb 



w;3cosh 2 [^] 
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Fig.(^]) shows /C(t) vs. t for temperatures T = 0,1,5,10. Again the oscillations are a conse- 
quence of the bound state contribution, and as in the Sine-Gordon case we find that despite the 
fact that in the high temperature limit the kernel becomes very localized in time, the total integral 
drKij) = preventing a representation of the noise-noise correlation function as a delta func- 
tion in time even in the high temperature limit. The "color" in the noise-noise correlation function 
is enhanced by the coupling to the continuum via the bound state which is also responsible for the 
strong oscillatory behavior of the real-time correlation function. 



5 HARD THERMAL LOOPS VS. CLASSICAL LIMIT 

The high temperature limit corresponds to T » m with m the meson mass. However we are 
restricted to the dilute limit in which the treatment of isolated domain walls is meaningful. Because 
the kink density is suppressed by an Arrhenius activation factor [H] 

iVfc~e-f, (112) 
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the study of the high temperature limit for the dynamics of isolated domain walls requires that the 
temperature range be such that 

TTi 

m « T « Mk— . (113) 
A 

For weak coupling A << 1 there is a wide temperature range in which the high temperature and the 
dilute kink gas approximation will be reliable. In order to understand the high temperature limit 
it is convenient to separate the loop integrals into the "soft part" in which both the integrated and 
transferred momenta are "soft", i.e. k,Q « T and the "hard" part, in which the loop momentum 
k is 0(T). Since the structure factor S(Q) is strongly suppressed for Q » m, the transferred 
momentum Q is always of order of m << T hence it is always "soft". The "hard" k region 
with "soft" transferred momentum is the domain of validity of the hard-thermal loop resummation 
programme of Braaten and Pisarskip8[|. 

Then for the soft region of the remaining k integral, we can replace the occupation factors 
n k ~ T/u) k . This soft region therefore gives the classical contribution to the kernels T,,T. A 
simple WKB analysis of the continuum solutions for both cases considered, reveals that the matrix 
elements D ktP fall off as ~ S(p + k)/k in the limit in which Q = k + p ~ m ; k — > oo. This simple 
analysis is confirmed by the exact expression for the matrix elements T k}P (see equations 8^JlQ5| ) 



which in this limit (Q = k + p ~ m, k — > oo) behave as S(Q)/k. Therefore in the meson loop, the 
matrix elements yield a contribution of 0(l/k 2 ) in the hard thermal loop limit, for which a simple 
scaling analysis reveals a large temperature behavior of 0{1/T). Hence we see that in the 1 + 1 
dimensional case the hard-thermal loop limit yields a subleading contribution as compared to the 
classical contribution from the soft region. This is a consequence of the small phase space available 
for the loop integrals in 1 + 1 space-time dimensions. 

This analysis allows us to conclude that the high temperature limit is dominated by the classical 
contribution with a lineal dependence on temperature in the regime T >> m. This behavior is 
clearly displayed in figure (Q) that shows the integrals I\ ; I2 with 

h = I dk 3 — ^-3 n k 



00 



k^k 



dk ,,3 ,3 n fc> ( 114 ) 

-00 ^Q-k^k 



corresponding to the contributions from Si ; S2 to the self-energy kernel at r = for Sine Gordon 
theory, with similar results for c/> 4 . We clearly see that for T > 2m the temperature dependence 
becomes lineal. Furthermore we have numerically checked that most of the contribution in this 
regime arises from the "soft" region of the loop momentum k < m and Q < m. 

Combined together the hard-thermal loop analysis and the numerical evidence lead us to con- 
clude unambiguously that the high temperature limit of the self-energy kernel is dominated by the 
classical finite temperature contribution. 

This analysis also holds for the noise-noise correlation function (since the same matrix elements 
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contribute to them). However for the noise- noise correlation function there is an extra factor of 
the Bose occupation factors in the integrals. This results in one extra power of temperature in 
the "soft" region while the temperature dependence from the hard-thermal loop region is mostly 
unaffected by the extra Bose factor. Therefore we conclude that the self-energy kernel is 0(m 2 T) 
and the noise-noise correlation function is 0{m 2 T 2 ) in the high temperature limit. This is in accord 
with the classical Fluctuation-Dissipation theorem in which the noise-noise correlation function has 
an extra power of temperature compared with the dissipative contribution. 

At long times the contribution from the two meson cut gives a rapidly oscillating contribution 
leading to a rapid fall off of the time dependence. On the other hand, the contribution from 
the Landau damping cut gives the leading contribution at long times because the discontinuity 
has support at very low frequencies and completely dominates long time behavior. Therefore 
we conclude that the long time, high temperature behavior in the dilute kink limit is completely 
dominated by classical finite temperature dynamics and dominated by the contribution from Landau 
damping. 

6 HIGHER ORDERS AND HIGHER DIMENSIONS 

At two-loops and higher orders we expect that collisions will provide a non vanishing static friction 
coefficient and result in an exponential relaxation of the velocity in some time regime. However, the 
contribution from these terms will be of higher order in coupling (m/M) and therefore there will be 
a competition between the time scales associated with lowest order relaxation via off-shell Landau 
damping and the higher order collisional relaxation leading to an exponential fall-off. Therefore we 
anticipate several different relaxational regimes with wide separation of the time scales for weak 
couplings and in the dilute regime. 

In 3+1 dimensions for degenerate scalar potentials the situation is clearly more complicated. 
The zero mode from translation invariance now gives rise to two-dimensional massless degree of 
freedom corresponding to small local distortions perpendicular to the (planar) wall. These are the 
capillary waves fluctuations of the interface that will dominate the long- wavelength small frequency 
dynamics. We expect to report on further studies of higher order collisional relaxation as well as 
new phenomena in 3+1 dimensions in the near future. 

7 CONCLUSIONS AND FURTHER QUESTIONS 

We have studied the non-equilibrium dynamics of domain walls in 1 + 1 dimensional scalar field 
theories at finite temperature in the dilute regime. We obtained the real time equations of motion 
for the expectation value of the collective coordinate and also the quantum Langevin equation to 
lowest order in the weak coupling (adiabatic) expansion. Two specific models were studied: (f> A and 
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Sine-Gordon scalar field theories providing detailed analytic and numerical studies of the equations 
of motion and a Markovian approximation to it. 



To lowest order in weak coupling we found that the real-time equation of motion involves a 
non-Markovian self-energy kernel and that the static friction coefficient vanishes. However, there 
is dynamical friction which is a result of the memory effects in the self-energy and is associated 
with two different types of two-meson processes: spontaneous and induced two-meson creation and 
annihilation and scattering off in medium mesons. The second type processes only occur at finite 
temperature and lead to Landau damping. 

We studied the Markovian approximation and shown numerically that this approximation is 
unreliable in a wide range of temperatures. 

The quantum Langevin equation was obtained by computing the influence functional obtained 
by tracing out the meson degrees of freedom to the same order in the adiabatic expansion. We 
found that the dissipative kernel and the noise correlation function obey a generalized form of 
fluctuation-dissipation relation but that a Markovian limit is not available, the noise is Gaussian, 
additive but colored. 

The high temperature limit in the dilute regime was studied in detail by analyzing the "soft" 
and "hard" contributions to the self energy and noise kernels. We find that the hard contribution is 
suppressed at high temperature because the matrix elements fall of as an inverse power of the hard 
momentum. The small one dimensional phase space leads to a suppression of the hard momenta 
and therefore the leading contribution at high temperature arises from the "soft" region with 
momenta of the order of the meson mass yielding the classical result in the high temperature limit. 
Furthermore the long time dynamics is completely determined by the Landau damping processes in 
the medium leading to the conclusion that to the order studied the long time dynamics is completely 
determined by classical Landau damping. 

We have restricted our study to a perturbative expansion which already showed the complicated 
nature of the problem even at lowest order. Pursuing a higher order calculation in perturbation 
theory will clearly be a major task. When the velocity of the soliton becomes very large one 
must abandon the approach advocated in this article and pursue a non-perturbative approach that 
accounts for strongly non-linear processes. We are currently implementing such an approach in 
terms of a self-consistent variational method [49] and expect to report on new results in the near 
future. 
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A REAL-TIME MESON CORRELATION FUNCTIONS 



In this appendix, we will calculate the Green's functions which are defined in eqn.(65) in terms of 
the vertex given by eqn. 



Applying Wick's theorem and eqn.(|3q), it is a matter of straightforward algebra to find the 
following results: 



p,k^0 I 



-i(oj p +u) k )(t— t 



') (n p n k + 9(t - t')(l + n p + n k )J 



, r) rp{A) rp{A) 

L1 pk 1 -p-k 



+ e+t^M (n p n k + 9(t' -t)(l+n p + n k ) 
i{u p -u k ){t-t') ( npHk + Hp Q{t , _ t)+nk 6{ t _ t >) 



G> (t, t')9(t - t') + G< (t, t')6{t' - t) 



G—(t,t>) 



p,k^0 I 



I O rp(A) rp(A) 

L± pk ± -p-k 



-i(u>p+u>k)(t—t 



') (n p n k + 9(t' -t)(l + n p + n k )) 



+ e+'^^'C- 1 ') (n p n k + 9{t - t'){l + n p + n k ) 
i(up-u k )(t-if) ^ npnk + nk e( jJ - t )+n p 9{t - t') 



G> (t, t')9{t' -t) + G < (t, t')9{t - t') 



G + -(t,t>) 



2 T,\ T S )T - 

p,k^0 I 



(S) rp(S) 

p—k 



e - i{ u p+Uk )(t-t') HpHk + e +i(^k)(t-t')^ npnk + np + nk + 1 



, r) rp(A) rp(A) 

L1 pk 1 -p-k 



<"p-"k)(t-t') np (i + nk ) 



-G<(t,t^ 
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e- i( ^+^ )(< - < ' ) (n p n fe +n p + n k + l) + e +i(^+^)(t-t') ^ 



I 9 rp\A) rp(A) 

L1 pk 1 -p-k 



-i(w p -w k )(t- 



n k (l + ra p ) 



G > (t,t') = -G < (t',t). 



(115) 



These Green's functions satisfy the following relation 

C++ + G + G+- + G~+ = (116) 

which is a consequence of unitary time evolution [42]. 

Furthermore, using the antisymmetry property of the matrix elements T^ k one finds that 

G+-(t,t') = (G-+(t,t')y (117) 

The Green's functions G ++ (t,t') ; G~~(t,t') can be written in terms of G^ (t,f) and its complex 
conjugate, therefore we see that there is only one independent Green's functions (and its complex 
conjugate). 

B CALCULATING K x (t - t') AND K(t - f) 



Performing the coordinate transformation in eqn.(pq), the influence- functional becomes 



T[x, R] 



exp 



1 



dtdt' 



R(t) R(f) 



G ++ {t,t') + G (t,t') - G + -(t,t') - G- + (t,t') 



+ ^R(t)x(t')(G ++ (t,t') - G—(t,t') + G + -(t,t') - G- + (t,t r 
+ h{t)R(t')(G ++ (t,t') - G—(t,t') - G + ~(t,t') + G" + {t,t') 



(118) 



Integrating the linear term in R by parts once and the quadratic term twice, the influence- 
functional can be cast in the following form 



F[x,R] =exp\ - I dtdt' R(t)K l {t-t')±(t') - R(t) K{t - t') R(t') 



(119) 



where 



1 d 

KAt-t') = ~ 
v ' 2dt 



(g++(M') - G—(t,t') + G + ~{t,t') - G- + (t,t' 
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+ (c ++ (t',t) - G (t',t) - G + ~(t',t) + G~ + (t',t 

i§- t [G > (t,t')-G<(t,t')]e(t-t>) 
i d 2 



G ++ (t,t') + G—(t,t') - G + -{t,t') - G- + {t,t') 



\^[G>(t,t') + G<(t,t') 



(120) 



(121) 



Substituting the values of the Green's functions from eqn. flllq ) in the above equations, one 
obtains the expressions for K\(t — t') and K (t — t') in eqns.(69) and (|70|). 



In the case that there are bound states other than the zero mode, such as the case of </> 4 the 
sum in eqn.(|32j) runs over all bound and scattering states, i.e. 

1 



F[a\a] 



2i 



dpdk.l—D p k 



,t J 



Ofcflp — a_ k a_ p + a_ k a p — a_ p ak 



1 



-Dbk 



OfcOfe - a^ k a\ + a)_ k a b - a\a k 



2i 



E / dk 



-Dkb 



dbdk - aW-k + a b a k - a -k a b 



a,b 



a b a a ~ 4°a + a fe a a ~ a l a fe 



(122) 



where the indices a and 6 stand for summation over discrete bound states and p and k stand for 
summation over continuum scattering states. The models which we considered in this paper have 
at most one bound state, that is the case in the cp theory. In this case, the last term will not 
contribute since D bb vanishes. Thus for only one bound state, eqn, (|122j ) can be written as 



F[a\a] 



dpdk [T^ {a p a k - al p ai k ) + T^ k } (ai p a k - a j _ k a p 



+ / dk 



T b f (a k a b - al k a\) + (a\ 



k a b — a b ak 



(123) 



where the matrices T^.? and for scattering states are given by eqn.(|33|) and if one of the states 



pk 

is a bound state, then 



AS) 

L bk 

n(A) 
L bk 



1 

2i 
1 

2? 



1 Mb 
^k 

F^b 



+ 



\Uk_ 

UJ b 
UJ b 



D b k 

Dbk ■ 



(124) 



In the Sine-Gordon theory, the last two terms in eqn.( |123| ) do not contribute since in this theory 
there are no bound states other than the zero mode and the Green's functions are given by eqn. (|ll5|) 
but with integration over p and k instead of the summation. 
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In the 4> case, to lowest adiabatic order the contributions from the bound and scattering states 
decouple. This implies that the Green's functions will have a contribution from the bound state 
which is given by the same expression as that of the scattering states, with p — > b, but multiplied 
by a factor of 1/2 since the bound state wave function is chosen to be real. 



C GENERALIZED FLUCTUATION-DISSIPATION RELATION 



The functions 



G>(t-t') = (F(t)F(t')) 
G<(t-t') = (F(t')F(t)), 



(125) 
(126) 



(where F is given by eqn.(|32|)) admit a spectral representation, and their Fourier transforms in 
time, g > (ui) ; g < (u) obey the KMS condition |50|| 

(127) 



g<{u) = e-^g>{u). 



From this expression we find that ki(tu), the Fourier transform in time of the kernel iKi(t — t') 
2£ m (i — t') is given by 



fci(w) = 2 



1 - e-P"' 



2tt 



uj — uj + te 



leading to the imaginary part 



Im\ki(ijj)} = — u>g > (u>) 



1 - e 



(128) 
(129) 



On the other hand the kernel that determines the noise- noise correlation function K(t — t') has 
a Fourier transform given by k(oj) with 

,2 



-%-l g > {u))+g <^)]=-\g>{u) 



1 + e 



-/3uj 



— Im[ki (uj)] coth 



(5u 



(130) 



The relation between the Fourier transform of the noise-noise correlation function and the imaginary 
part of the self-energy is the generalized Fluctuation-Dissipation relation ^ ]. 
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Figure Captions 



Figure 1 The non-equilibrium one-loop contributions to the self energy. The upper two con- 
tributions correspond to emission-annihilation of two mesons. The lower two correspond scattering 
off in-mcdium mesons and responsible for Landau damping. 

Figure 2 Contour in the complex s-plane for the inverse Laplace transform. 

Figure 3 The functions T(t) and S(r) for temperatures T = 0, 1.0, 5.0 and 10.0 for Sine- 
Gordon theory. 

Figure 4 Numerical evaluation of the velocity of the kink for j = ; vq = I for temperatures 
T = 0, 1.0, 5.0 and 10.0 in Sine-Gordon theory. 

Figure 5 The correlation function /C(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the 
Sine-Gordon theory. 

Figure 6 The functions r(r) and S(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the 
(/> 4 -theory. Contributions from bound and scattering states are displayed separately. 

Figure 7 Numerical evaluation of the velocity of the kink for j = ; vq = 1 for temperatures 
T = 0, 1.0, 5.0 and 10.0 in cf) 4 theory. 

Figure 8 The correlation function /C(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the 4> A 
theory. 

Figure 9 Integrals I\ ; I2 corresponding to the contributions from Si ; £2 to the self-energy 
kernel at r = vs. T in the Sine Gordon theory. 
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